Effect of disorder and electron-phonon interaction on interlayer tunnelling current in 

quantum Hall bilayer. 
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We study the transport properties of the quantum Hall bilayers systems looking closely at the 
effect that disorder and electron-phonon interaction have on the interlayer tunnelling current in the 
presence of an in-plane magnetic field Bu . We find that it is important to take into account the 
effect of disorder and electron-phonon interactions in order to predict a finite current at a finite 
voltage when an in-plane magnetic field is present. We find a broadened resonant feature in the 
tunnelling current as a function of bias voltage, in qualitative agreement with experiments. We 
also find the broadening due to electron-phonon coupling has a non-monotonic dependence on B\\ , 
related to the geometry of the double quantum well. We also compare this with the broadening 
effect due to spatial fluctuations of the tunnelling amplitude. We conclude that such static disorder 
provides only very weak broadening of the resonant feature in the experimental range. 



I. INTRODUCTION 

Over the last fifteen years quantum Hall bilayer sys- 
tems (QHB) have been extensively studied since they are 
one of the few systems that show macroscopic evidence 
of quantum coherence. The richness of the physics of the 
QHB has attracted the attention of both theoretical^— 
and experimental studies^—. This has led to rather 
rapid progress in the area. The bilayer consists of two 
parallel two-dimensional electron layers in a double quan- 
tum well closely separated by a distance d and subjected 
to a magnetic field perpendicular to the plane of the 
layers B±. In this paper, we will focus on the case 
when each layer is a half-filled Landau level: filling factor 
Vy = 1/2 = 1/2. If the separation between layers is large 
they behave as two independent Fermi liquids and no 
quantum Hall effect is observed. When the distance be- 
tween the layers becomes comparable with the magnetic 
length (d ~ £b) the system undergoes a phase transition 
from a compressible state at large d 3> £b to an incom- 
pressible state small d <C Ib- In the incompressible state, 
the system as a whole exhibits the v = 1 quantum Hall 
effect even when interlayer tunnelling is negligible. This 
transition is driven by Coulomb interactions between the 
layers. The ground state in the quantum Hall regime is 
believed^ to have a broken U(l) symmetry which leads 
to spontaneous interlayer phase coherence. This ground 
state can be described as a pseudospin ferromagnet£ us- 
ing a pseudospin picture equating electrons in the upper 
(lower) layer with pseudospin "up" ( "down" or as 
an excitonic superfluid 4 ^. 

A series of remarkable experiments have probed the 
existence of this coherent phase. They show evidence for 
interlayer coherence with a linearly dispersing Goldstone 
mode 11 and counterflow superfluidity 9,14 ' 15 and drag Hall 
voltage 16 . One piece of the experimental evidence of in- 
terlayer coherence is a sharp peak in the tunnelling cur- 
rent for small bias (between 10 and 100/xV) and low tem- 
perature. Another characteristic of the QHB that indi- 
cates the existence of phase coherence between layers is 



the sensitivity of the QHB to the presence of an in-plane 
magnetic field. The sharp peak in the tunnelling cur- 
rent at small bias is suppressed when a magnetic field is 
applied parallel to the plane of the layers. At the same 
time, experiments show a 'dispersive' feature in the tun- 
nelling current at a voltage which evolves linearly with 
the in-plane field Bu . This has been interpreted as the 
excitation of the Goldstone mode of the excitonic super- 
fluid at wavevector Q and energy eV given by 



Q = 



2-KB^d 



eV = HvQ 



(1) 



where v is the velocity of the collective mode and 4>o — 
hc/e is the flux quantum^ 1 -. 

The tunnelling current can be computed using as a 
perturbation the interlayer hopping matrix element tj_. 
For a homogeneous system in the absence of an in-plane 
field, Jack et ai 6 - found a tunnelling current proportional 
to 1/V at zero temperature, consistent with the observa- 
tion of a region of negative differential tunnelling conduc- 
tance at low temperature. However, the same calculation 
gives a delta function at eV — vQ in the tunnelling cur- 
rent in the presence of an in-plane field. Although the 
position of this feature is consistent with experiments, 
the experimental data does not exhibit a sharp feature 
even at the lowest temperatures. In this paper, we will 
explore sources for a finite linewidth of this feature. We 
find that two mechanisms should be dominant at low 
temperatures: tunnelling disorder and electron-phonon 
coupling. We will not discuss the role of vortices which 
could be nucleated at zero temperature by strong charge 
disorder or by thermal activation. 

This paper is organised as follows. In the next section, 
we review the methodology for computing the tunnelling 
current as originally used by Jack et al^. In Section IlIIl 
we discuss the effect of electron-phonon interactions on 
the tunnelling current. We will calculate how the magni- 
tude and the width of the dispersive feature in the tun- 
nelling current is affected by the electron-phonon cou- 
pling. In Section IIV1 we investigate the effect that an 



2 



inhomogenous tunnelling amplitude over the sample has 
on the tunnelling current. We present the conclusions of 
the paper in Section Ivl 



II. THEORETICAL DISCUSSION 

We will adopt the pseudospin picture of the quantum 
Hall bilayer, labelling single-particle states in the up- 
per layer as S z = +1/2 and states in the lower layer 
as S z = —1/2. We will review this framework in this 
section. In this picture, the system is a pseudospin ferro- 
magnet with an easy-plane anisotropy. Furthermore, we 
will work with the large-S* generalisation of this model on 
a lattice which corresponds to coarse-graining the system 
by treating ferromagnetic patches containing S electrons 
as lattice sites containing a large spin S. The Hamilto- 
nian can be written as: 



Ho = ~J2( S t S 7 + h - c -) + T 5>* Z ) 



D 



(2) 



where Sf and S~ are the pseudospin raising and lowering 
operators on site i of a square lattice, J is the exchange 
interaction, D represents a local capacitative energy for 
charge imbalance. In the presence of tunnelling across 
the bilayer, the Hamiltonian becomes: 



H = H n - A 



(3) 



where Aj^ is the inter layer tunnelling matrix element in 
the absence of in-plane field and Xi is the ^-coordinate of 
the spin i. We have chosen the gauge such that, in the 
presence of an in-plane magnetic field B\\ , the tunnelling 
matrix acquires a phase that varies spatially in the re- 
direction with periodicity 2it/Q with Q as defined in ([1}. 

The large-S* treatment of this model corresponds to 
taking S to a large value while keeping JS and DS con- 
stant so that the three energy scales for exchange, inter- 
action and tunnelling scale in the same way with S. We 
note that DS > JS 3> Aj_ in the bilayer system. 

The pseudospin can be written in terms of phase and 
S z operators as 



(S-S z ) 1/2 e-^(S + S z ) 1/2 



(4) 



Semiclassically, 4>i gives the azimuthal angle of the pseu- 
dospin projected onto the a;y-plane in spin space. They 
are canonical conjugate variables: [S£,<frj] = iSij. We 
can take the continuum limit and integrate out the S z 
(charge imbalance) fluctuations to arrive at a phase-only 
actional: 



S = 2Sp s d 2 rdt 



2u z 



)2_i|V^| a + icosfa-Qa;) 



(5) 

where v 2 = DSJSi 2 B /Ti 2 , p s = JS is the spin stiffness 
and £ = (47T/9 s /Ax) 1 ^ 2 ^b is the Josephson length which 



gives the length scale over which counterflow currents 
decay due to tunnelling across the bilayer. If we further 
include a bias of V across the bilayer, the only change to 
the action is that the cosine term in the action becomes^ 
cos(4> — Qx — eVt/h). In this model, all the dynamics of 
the system depends on the phase <fi. Spatial gradients in 
the phase correspond to counterflow in the two layers. 

In the absence of tunnelling, this Lagrangian repre- 
sents an easy-plane ferromagnet with spontaneously bro- 
ken symmetry, i.e. a spatially uniform phase. In the pres- 
ence of an in-plane magnetic field, the ground state de- 
velops spatial variations in the phase field 1 -' 1 -. At small 
Q, the phase field increases linearly with the Aharonov- 
Bohm phase: ~ Qx. At large Q, the phase field cannot 
follow the Aharonov-Bohm phase and has only small os- 
cillations: ~ sin(Qx)/Q 2 ^ 2 where £ is the Josephson 
length. The transition occurs at Q£ ~ 0(1). The Joseph- 
son length is estimated to be of the order of microns. 
Experimental values 11 of the in-plane field give Q£ > 10. 
Therefore, the spatial fluctuations of the ground-state 
phase field are small. In our following calculation, we will 
consider quantum fluctuations as perturbations around 
the uniform = state. 

Our calculation of the quantum fluctuations in the sys- 
tem starts by the system in the absence of tunnelling. 
Then we will introduce interlayer tunnelling in perturba- 
tion theory. In the absence of tunnelling, the Aj_ = 
Hamiltonian ([2]) can be diagonalized 



H 



£ a = 



[JS<t>l%(DS + ?l%)]W 



at = \ / - 



S z 



(6) 



where the coherence factors are (u q + Vq) 2 = (w q — 
v^y 2 = (1 + DS/JSq 2 £ 2 B )^ 2 . This means that the ele- 
mentary excitations in the pseudospin lattice system, cre- 
ated by the operator a q are long-lived pseudospin waves 
with energy £ q at wavevector q. At long wavelengths, 
the pseudospin waves have a linear dispersion, e q ~ hvq, 
with pseudospin wave velocity v - (JSDS)^ 2 l B /h. 

The operators S z , e ±l< ^ and 5* all involve the creation 
and annihilation of pseudospin waves. Most significantly, 
interlayer tunnelling, S ± , causes decay of the pseudospin 
waves at all wavelengths. This can be seen by examining 
the S^~ operator that represents electron tunnelling in 
the pseudospin language. From equation (j4j, we see that 
it creates perturbations in the phase 4>i and Sf . This in- 
volves the creation and annihilation of pseudospin waves 
(see equation ([6])). This is a consequence of the fact that 
tunnelling breaks the global U(l) phase invariance of the 
system so that Goldstone's theorem no longer protects 
the long-wavelength pseudospin waves from decay. 

Let us consider now the interlayer tunnelling current 
in the presence of an interlayer bias V. The tunnelling 
current at site i of the lattice is given by the operator 
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iA ± (e~ l ^- eVt /^S- -e l(Qx >- eVt / n S+)/2h. Therefore, 
in the continuum limit, the expectation value of the in- 
terlayer tunnelling current ia^: 

w-^/^W'-o-tt)) <7) 

where the expectation value is taken with respect to the 
full Hamiltonian H as defined by <j3j> - In a perturbativc 
treatment of the interlayer tunnelling in the Hamiltonian 
H , we can treat the tunnelling term in first-order pertur- 
bation theory. This gives^ a dc current proportional to 
Ai: 

(8) 

where Q = (Q,0), C(r,r) = ((Te^^e"^ ' )^ - h.c.) 
evaluated at zero tunnelling with Hamiltonian Hq. We 
can see that we are calculating the response of the system 
at wavevector Q and frequency eV/%. It can be shown 
that 

ln(Te ^) e ^(r',o) )Q = ^/|^|J q ^(r,r',r)| 2 G5 

(9) 

where J q , w (r, r', r) = e^"--^) - e l * r ' and G^l = 
(uj 2 /v 2 — q 2 )^ 1 is the propagator for the phase fluctu- 
ations at zero tunnelling. 

Jack et ai* analysed the tunnelling current at zero in- 
plane field. They showed that the perturbative calcula- 
tion above can be understood in terms of the generation 
of finite-momentum pseudospin waves via the decay of 
the q = mode. Technically, this is an interpretation 
of {8j as a Taylor expansion of the exponential in the 
definition of C(r,t). Each term in the Taylor expansion 
involving In 0-fields represents the generation of n pseu- 
dospin waves. In the absence of an in-plane field, there is 
no decay of the q = mode to a single pseudospin wave 
with non-zero wavevector because of momentum conser- 
vation. The most important decay channel is then the 
generation of a pair of pseudospin waves with equal and 
opposite momenta. However, in the presence of an in- 
plane magnetic field Bh, the vector potential provides a 
momentum of Q, as given by equation ([1]), to the pseu- 
dospin system, as can be seen in the Hamiltonian ([3]). 
The generation of a single pseudospin wave at wavevec- 
tor Q is now possible and this is the leading contribution 
in orders of l/S. (In the perturbative formulation, this 
can be seen mathematically in <j8j) which tells us that we 
need to calculate the Fourier component of C at wavevec- 
tor Q.) For a homogeneous system, this calculation gives 
a current which is a delta function at eV/Ti = vQ: 

In this work, we investigate possible sources of line broad- 
ening for this peak at non-zero in-plane magnetic field. 



We will focus on effects which do not vanish at zero tem- 
perature. In order to obtain a finite linewidth, we find 
that the pseudospin waves need elastic or inelastic scat- 
tering. We will discuss elastic scattering due to disorder 
in section HVl In section Hill we will study inelastic scat- 
tering. This can arise from the generation of photons 
or phonons. The two mechanisms are similar. However, 
we will see that the energy of photons involved in this 
process will be much higher than the energy of the pseu- 
dospin waves. This means that the process can only be 
virtual. It can at most alter the dispersion relation of 
the pseudospin waves but cannot cause decay. There- 
fore, we will focus on the electron-phonon interaction in 
this work. 



III. PHONON GENERATION 

In this section, we introduce interactions between 
phonons and pseudospin waves and study how the tun- 
nelling current between layers is affected by the introduc- 
tion of these interactions. 

When electrons tunnel across the bilayer, the electron 
density changes, perturbing the core ions on the AlGaAs 
of the tunnelling barrier and thus creates phonons in the 
three-dimensional system in which the quantum well is 
embedded. For the range of values of the bias voltage 
used in the experiments, the most important interaction 
between ions in the host material and the tunnelling elec- 
trons is the deformation potential interactio n 20 ^ 1 with 
the acoustic phonons. We neglect optical phonons be- 
cause they are at energies high compared to the electron 
energy at the experimental range of bias voltage. 

The phonon Hamiltonian is given by: 

-ffphonon = ftw k ^a^ak + (11) 

where k = (q, k z ) is a three-dimensional wavevector, q = 
{q x ,qy)i is the phonon frequency spectrum and Ok and 
are the phonon annihilation and creation operators. 
We are discussing physics at wavelengths long compared 
to the lattice spacing of the substrate. So, we will use a 
simple linear dispersion for the acoustic phonon: = 
s|k|. The electron-ion interaction takes the form: 

/ * \ 1/2 
He-ion = D ci 5 |k|p(k)(a k + a T _ k ) 

(12) 

where p lorl is the ion mass density, v is the volume of the 
three-dimensional solid and p(h) is the three-dimensional 
electron density. 

The electron density perturbation caused by tunnelling 
involves charge imbalance across the bilayer. It is con- 
venient to express the density perturbation at position 
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R 



(r, z) in terms of the z-component of the pseudospin: 



p(R) 



4n£% 



Pt( z )+Pl( z ) + 



S z (r) 



(j>f(z) -pi(z)) 



(13) 

where p(Et) is the three-dimensional electron density and 
Pf(i)( z ) i s the one-dimensional density profile in the z- 
direction for electrons in the upper (lower) layer, normal- 
ized to J p^m(z)dz — 1. We have scaled the density by 
S in the spirit of the coarse-graining idea of the large-.!? 
generalisation of this model, as discussed at the start of 
section HH It is easy to check that the above expression 
produce the expected three-dimensional density profile 
for the cases of to a balanced bilayer S z = and com- 
pletely imbalanced bilayer S z — ±S. 

For simplicity, we approximate the electron wavefunc- 
tion ip(z) in the lowest subband of the single quantum 
well centered at z = by a Gaussian wavepacket: 



r/>(z) 



rV4 



(14) 



where w is the width of the well. (A factor of 7r 2 has been 
inserted into the exponential so that the form matches 
the amplitude and width of the actual wavefunction in 
the square well.) Using this approximation, the electron 
density in the Fourier representation is given by 



P(k) 



k z d 
2 



-k 2 w 2 /47T 2 QjZ _|_ 



(15) 



where d is the separation of the centres of the two quan- 
tum wells in the z-direction. We have dropped terms 
that do not involve the pseudospin degrees of freedom, 
i.e. terms that do not change as a result of a tunnelling 
event that creates a charge imbalance across the bilayer. 
Note that the part proportional to S z vanishes as k z — > 
because this part represents charge imbalance and does 
not contribute to a uniform change in density in the z- 
direction. In fact, the form of the fc z -dependence in p(k) 
reflects the geometry of the double quantum well. There 
are two contributions to the geometrical form factors in 
this expression, analogous to the Fraunhofer diffraction 
pattern from a double slit in optics. The sinusoidal de- 
pendence on k z d/2 comes from the convolution of the two 
density profiles centered at ±d/2 and the Gaussian comes 
from the electron density profile in the z-direction. From 
equations (@}, © and (flip)) , we see that the electron- 
phonon interaction in equation (I12[) includes the decay of 
one pseudospin wave with momentum q into one phonon 
with momentum (q, k z ) (or vice versa). This is depicted 
in Figure [TJ There is no conservation of momentum in 
the z-direction because translational symmetry is broken 
by the presence of the double quantum well. This decay 
of the pseudospin wave into a phonon can only occur as 
a real process if energy is conserved: 



vq 



;vV + kl 



(16) 



FIG. 1: Diagrammatic representation of the decay of one 
pseudospin wave with momentum q into one phonon with 
momentum q, k z . 



This requires the pseudospin wave velocity to be higher 
than the phonon velocity: v > s. This is the case we con- 
sider here. We estimate v ~ 1.4 x 10 4 ms _1 from the data 
of Spielman et al^- and s ~ Q.lv for the sound velocity in 
the heterostructure^i. On the other hand, if the phonon 
velocity s was greater than the pseudospin wave veloc- 
ity v, then the process where a phonon of greater speed 
(or indeed a photon) can be emitted and re-absorbed is 
only a virtual process. This will alter the pseudospin 
wave spectrum but does not contribute to the tunnelling 
current. 

We will now investigate more quantitatively how the 
pseudospin-phonon interaction affect the system. We 
proceed by integrating out the phonons to obtain an effec- 
tive action. This can be performed because the electron- 
phonon coupling is linear in the phonon field. We obtain 
a retarded interaction for S z , the charge imbalance across 
the layers: 



2Sv* fd^qdu 2 
6cp ~ 8n£%pJ (2tt)3 ADll ^ b ^\ 



where 



n q , w = / dk. 



sin 2 (/M/2)(q 2 + kl)e ~ 2k " w2 / A ^ 



(17) 



(18) 



and Xrj = (D c f/~h) 2 (p s / p ion Sv 2 ). This adds to the on- 
site instantaneous interaction D. Integrating out S z , we 
obtain a modified phase action at zero tunnelling: 



So 



(2tt) 3 7 V 1 + Aijll q , u 



\K^\ 2 (19) 



The decay of the pseudospin wave into a phonon is 
encoded in the imaginary part of n q . w . It is non-zero 
when uj > sq: 



i m n r 



co 2 sm 2 (d^oj 2 - ^V/2s)e" 2w > ~ s « ^ s 



2 9 

gig* 



( 20 ) 

For u) < sq, the imaginary part of ITq.^ is zero. The real 
part of Ilq w w iH shift the position of the peak of the 
tunnelling current with respect to the voltage. 

The leading term on the expansion of © in powers of 
1/S corresponds to the decay of one q = mode into one 
finite-momentum spin wave at wavevector Q = (Q,0). 
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a) 




q=0 



FIG. 2: Diagrammatic representation of: a) the decay of one 
q — pseudospin wave into a pseudospin wave with finite 
momentum Q. b) the decay of one pseudospin wave with 
momentum q — into a finite- Q pseudospin wave which in 
turn decays into a phonon. 



The contribution to the current from this term is 

2 



I = Se 



( A ± L 

\4tt?i£ 2 b 



Re/ dr du;e teVT/h cos(wr) x 

JO J -oo 



ih 



with 



V/^-Q 2 ]. 



(21) 



(22) 



The main difference between the present Green's func- 
tion and the one used by Jack et al and Balents and 
Radzihovsky is that this decay is affected by the phonon- 
electron interaction (Figure Eb) while the process consid- 
ered by them is not (Figure [2^l) . This will give a broad- 
ened peak on the current centered at eV = vQ. Phys- 
ically, this is a consequence of the fact that, once the 
q — Q spin wave is coupled to the phonons, it is no longer 
a sharp resonance. Mathematically, we have to inspect 
Gq jU which includes the decay of one zero-momentum 
pseudospin wave with energy hu) into a phonon with the 
same in-plane momentum. The poles of Gq, w establish 
the relation between the momentum Q and the frequency 
uj and the integral over time in equation (|21[) gives the 
restriction for the pseudospin wave energy fuv to eV . The 
final expression for the current is then 



/ = -e 



\4irM 2 B 



hv 2 t(l + A p IlQ,n) 
2 Ps e 2il(n - eV/H) '■ 



(23) 



where Q is the pole of Gqq with a positive real part. 

A measure of the importance of the electron-phonon 
coupling is the dimensionless quantity (see Fig. [3]): 



Xjjlmll 



Q,vQ 



X D vQ . 2 
— o — sin 



Qdv 
~2s~ 



exp 



Q 2 w 2 v 2 
' 2ir 2 s 2 



This is formally small in the large-S* limit. Indeed, we 
can check that, even if we set S to unity, Xd^q.vQ re- 
mains small. Neglecting the geometrical form factors, the 
coupling strength is given by XdvQ/s 3 ~ 10 _1 for Q be- 
tween 10 and 20 /^m _1 (as probed in the experimentsii) . 
The geometrical form factors reduce this further. From 
Fig. 21 we see that the coupling is most appreciable 
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FIG. 3: Dimensionless coupling AoImflQ^Q as a function of 
Q = 27rB||d/0 o . (d =28nm, w =18nm, v/'s = 1 0, S = 1.) 



for wavevectors Q in the experimental ranged of 10 to 
20/im _1 (with d — 28nm and w — 18nm). This corre- 
sponds to Qdv /2s ~ n/2. On the other hand, the cou- 
pling can also vanish when Qdv /2s is an integer multiple 
of 7r. At these wavevectors, the electron and ion density 
oscillations are symmetric in the two layers and so are 
not excited by a tunnelling event which causes charge 
imbalance across the layers. 

The weak electron-phonon coupling means that we can 
expect the pole of the Green's function to be close to 
the original pseudospin wave energy: Q ~ vQ. We can 
approximate Uq,q by Hq, v q and 



n ~ vQ (i + \ D n QtVQ /2) 



(25) 



This corresponds to a spin wave decay rate of Tq = 
ADV<2|lmnQ,i,Q|/2. These results give us predictions for 
the height and width of the feature in the tunnelling cur- 
rent that disperses with the in-plane field Bu cx Q as 
reported by Spielman et aU^. From equation (|2"3")l . the 
maximum value of the current occurs at eV/Ti = KeQ ~ 
vQ. At eV/h = vQ, 



1 



1 



Ahp s Q 2 \ D lmRQ tV Q 



(26) 



The width AV of the peak in the current as a function 
of the bias V is similarly controlled by Imfi: 



eAV/h = vQX D \Irnn QtVQ \/2. 



(27) 



First of all, the area under the peak in the IV curve 
can be approximated by 7 max Al^. This is proportional 
to 1/Q so that this peak weakens as we increase the in- 
plane field. However, the evolution of the shape of the IV 
curve is a non-monotonic function of Q (see Fig. 2]). The 
peak is sharp when electron-phonon coupling is weak, 
i.e. away from Q ~ 10/im -1 . The peak is broad in the 
experimental regime because of the appreciable electron- 
phonon coupling as we noted above. Our results are 



6 



150 - 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

eV/ftClO 11 *- 1 ) 

FIG. 4: Tunnelling current as a function of the voltage eV/h 
for different values of Q — 2nB^d/cj>o (taken every I/an -1 ). 
From left to right, the first curve corresponds to Q = 10/zm - , 
and the last curve corresponds to Q = 20^im~ 1 . We use 
D ei = 9eV. 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 

eV/^ao'V 1 ) 

FIG. 5: Tunnelling density conductance dJ/dV as a function 
of the voltage eV/Ti for Q — lO/im -1 (left) and Q — 14/im -1 
(right). 



qualitatively consistent with the dispersive feature in the 
IV characteristic observed by Spielman et al 11 . This 
feature has been identified as due to the excitation of 
coherent excitations of the interlayer-coherent phase be- 
cause its position moves linearly with the in-plane mag- 
netic field. The width of the feature appears to have non- 
monotonic dependence on the in-plane field. It would be 
interesting to explore this dependence for a larger range 
of wavevectors in order to elucidate relaxation mecha- 
nisms in the system. 

We note that our theoretical results differ in absolute 
magnitude from the experiments. Our current values are 
nearly an order of magnitude higher than the experimen- 
tal values^ 1 -. This may be due to a strongly renormalised 
tunnelling amplitude at low energies due to fluctuations. 
We also predict sharper peaks than seen in experiments. 



The broadest peak we obtained (which does correspond 
to the range of in-plane fields studied experimentally), 
the theory gives a width AV of the order of 10% of the 
position of the peak V — hvQ/e. See Figure[5]for a direct 
comparison with the differential conductance. The dis- 
crepancy may be due to other relaxation processes which 
contribute additively to the broadening of the peak. (We 
will discuss the contribution from elastic scattering in 
the next section.) It could also be caused by a suppres- 
sion of tunnelling due to thermal 4 and quantum 6 phase 
fluctuations. 

In summary, we have found that electron-phonon scat- 
tering may have a measurable contribution on the broad- 
ening of the coherent feature in the tunnelling IV char- 
acteric in the quantum Hall bilayer. In our theory, this 
broadening arises from the finite lifetime of the collective 
density excitation of the bilayer at wavevector Q due to 
the decay of electron density oscillations into lattice vi- 
brations with the same in-plane wavevector. The effect is 
in fact strongest in the range of in-plane magnetic fields 
studied experimentally, giving a width AV of the order 
of 10% of the position of the peak V = hvQ/e. 

IV. EFFECT OF DISORDER 

We now consider the effect of disorder on the tun- 
nelling current. Charge inhomogeneity can be present 
in the system in the case of strong disorder from the 
random distribution of dopants. This will nucleate tex- 
tures in the pseudospins (merons), producing a ground 
state with random vorticity in the phase of counter- 
flow superfluid. This has been investigated by Eastham, 
Cooper and Lee22 who found suppressed tunnelling in 
the sense that the length scale for counterfiow current to 
leak across the bilayer by tunnelling becomes enhanced 
by at least an order of magnitude. 

In this paper, we will focus on a system of weak dis- 
order where the ground state of the system is free of 
vortices. This applies to relatively clean samples which 
have small charge variations even on the length scale of 
a few magnetic lengths. We will consider in particular 
spatial variations in the tunnelling amplitude across the 
bilayer as a consequence of not having perfectly flat lay- 
ers. In the pseudospin model, inhomogeneous tunnelling 
corresponds to inhomogeneous Zeeman coupling to the 
in-plane component of the pseudospin so that the Hamil- 
tonian is: 

H a =H -A ± J2 + <ny Qxi S? + h-c] (28) 

i 

where Oi represents the fractional fluctuation around the 
mean tunnelling matrix element Aj_. We will model the 
on-site fluctuation in <ii as Gaussian distribution with 
zero mean and variance a 2 . We will consider disorder 
which is short-ranged on the scale of the magnetic length. 
We can perform a perturbation theory around the zero- 
tunnelling ground state. After averaging over disorder, 
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the lowest order in perturbation theory (~ Aj_) would 
give a dc tunnelling current: 



I 



A ± L 

2hi% 



/poo 
d 2 R / dr C(R) x 

(l+^o-)e- l(Q - R+eVT/?i) (29) 



This is clearly not a physical result since this expression 
still contains the current Iq for the uniform system (see 
equation (ITUl) ) which give a delta- function response at 
eV/h — vQ. Therefore, we need to consider further cor- 
rections due to the disordered tunnelling strength. As 
we have learnt from the previous section, the scattering 
of the pseudospin wave at wavevector Q is important to 
consider. In the previous section, we considered scatter- 
ing due to electron-phonon interactions. Here, we need 
to consider scattering by the disorder. This process in- 
volves terms of the form 7 q ,Qa q o:Q in the expansion of 



the Hamiltonian (|28l) in orders of S. The matrix elements 
are of order S° and have the form: 



7q,Q = -2A_LCT q (u q + Vq)(u Q + Vq) 



(30) 



The decay rate of a pseudospin wave at Q can be calcu- 
lated from Fermi's golden rule as {2-k/K) J2q l7q,Q| 2 ^( e q — 
cq). We see that we require |q| = Q since these collisions 
are elastic. We can approximate 4(u q + u q ) 2 (u q + v q ) 2 by 
(DS/hvq)(DS/hvQ). Averaging over disorder, we obtain 
a decay rate of: 



A 2 ± a 2 (DSl B ) 



v 2 h 3 



HvQ 



A\a 2 
hp s QlB 



(3D 

Ps J 



for spatially uncorrelated disorder: aiUJ = a 2 Sij. 

Analogous to our previous treatment for electron- 
phonon coupling with a finite decay rate for the spin 
waves (|23|) . this scattering rate due to disorder gives a 
peak for the IV curve of width AV = KTq/e and height: 



Im ax I 



f A\ L 



\4n£ 2 B 



ev 1 L ev f p s \ 
8h Ps QTQ ~ ^128a 2 £ B \Ds) 



1/2 



This gives a field-independent peak current but a peak 
width that decreases with increasing Q. 

For disorder with a correlation length of C with the 
correlation function WiaJ = a 2 exp[— |r.; — we ob- 

tain: 



Tn = 



A\a 2 (DS_ 



TiPsQIb V Ps 



1/2 



c 2 



(1 + Q 2 C 2 ) 3/2 ' 



(33) 



As a function of the in-plane field, the peak width AV 
decreases and the peak current increases as we increase 
the in-plane field. This is observable if the in-plane field 
is large enough that Q( becomes large compared to unity. 

We can compare the width AV with the position 
V = HvQ/e of this resonant feature. For uncorrelated 
disorder, 



AV 



aA± 

PsQIb 



(34) 



In the experimental range, Q£b ~ 0.1, A±/p s ~ 10~ 3 
and we expect a < 1 for the validity of the perturbation 
theory. We see that this broadening is very weak. Conse- 
quently, the peak current is apparently very large as can 
be seen in our expression for the peak current (which is 
independent of Ax). Therefore, we conclude that spatial 
fluctuations in the tunnelling amplitude does not give rise 
to a strong broadening of the resonant IV peak. Con- 
versely, our results indicate that, since Tq ~ 1/Q, the 
broadening effect of this source of disorder is only ob- 
servable at much smaller values of the in-plane field. 



V. SUMMARY 

We have studied in this paper extrinsic sources of scat- 
tering for the collective excitations ('spin waves') of the 
counterflow superfluid, with the aim of understanding 
the broadened peak that disperses with in-plane field £?|| 
in tunnelling experiments. The IV peak then reflects 
the spectral weight of the spin waves at the wavevector 
Q = 2TrB ll d/<f>o. 

We have concentrated on disorder and electron-phonon 
interaction as sources of spin wave decay and hence a 
broadening of the IV feature. Interestingly, we found 
that the broadening is non-monotonic and is strongest in 
the experimental range of in-plane fields (Q ~ 10 /im -1 ) 
because the acoustic phonons emitted after the tunnelling 
event involve vibrations commensurate with the spacing 
between the two quantum wells. It will be interesting to 
investigate the linewidths more systematically to see if 
this monotonic evolution of the lineshape can be observed 
in experiments. 

Disorder also provides broadening. We have consid- 
ered fluctuations in the tunnelling amplitude across the 
bilayer. Our theory suggests that this provides only 
very weak broadening. Stronger disorder would nucle- 
ate charged quasiparticles in the system and is beyond 
the scope of this paper. This is discussed recently by 
Eastham et a&. 

However, the linewidths predicted are still nearly a fac- 
tor often smaller than the experimental results. As men- 
tioned above, this may be due to phase disorder due to 
thermal or quantum fluctuations. It will be very useful to 
have experimental measurements for a wider range of the 
magnetic fields and at lower temperatures to investigate 
whether electron-phonon or weak disorder are important 
in determining the lineshapes in this system. 



Acknowledgments 

We gratefully acknowledge the financial support of EP- 
SRC grant EPSRC-GB EP/C546814/01. 



8 



1 S.M. Girvin, A.H. MacDonald and J. P. Eisenstein, in Per- 
spectives in Quantum Hall Effects, ed. S. Das Sarma and 
A. Pinczuk (Wiley, New York, 1997). 

2 K. Moon, H. Mori, K. Yang, S.M. Girvin, A.H. MacDon- 
ald, L. Zheng, D. Yoshioka, and S.C. Zhang, Phys. Rev. B 
51, 5138 (1995). 

3 X.G. Wen and A. Zee, Phys. Rev. B 47 2265 (1993). 

4 L. Balents and L. Radzihovsky, Phys. Rev. Lett. 86, 1825 
(2001). 

5 A. Stern, S.M. Girvin, A.H. MacDonald, and N. Ma, Phys. 
Rev. Lett. 86, 1829(2001); M.M. Fogler and F. Wilczek, 
Phys. Rev. Lett. 86, 1833 (2001); Y.N. Joglekar and A.H. 
MacDonald,Phys. Rev. Lett. 87, 196802 (2001) 

6 R.L. Jack, D.K.K. Lee and N.R. Cooper Phys. Rev. Lett. 
93, 126803 (2004); Phys. Rev. B 71, 085310 (2005). 

7 J.P Eisenstein, G.S. Boebinger, L.N. Pfeiffer, K.W. West 
and S. He, Phys. Rev. Lett. 68,1383 (1992); J.P. Eisen- 
stein, L.N. Pfeiffer, and K.W. West, Phys. Rev. Lett. 69, 
3804 (1992). 

8 S.Q. Murphy, J.P. Eisenstein, G.S. Boebinger,L.N. Pfeiffer, 
and K.W. West, Phys. Rev. Lett. 72, 728 (1994). 

9 M. Kellogg, J.P. Eisenstein, L.N.Pfeiffer, and K.W. West, 
Phys. Rev. Lett. 90, 246801 (2003). 

LB. Spielman, J.P. Eisenstein, L.N. Pfeiffer, and K.W. 
West, Phys. Rev. Lett. 84, 5808 (2000). 

1 LB. Spielman, J.P. Eisenstein, L.N. Pfeiffer, and K.W. 
West, Phys. Rev. Lett. 87, 036803 (2001). 



S.M. Girvin, in proceedings of Les Houches summer school 
on Topological Aspects of Low Dimensional Systems, ed. A. 
Comtet, T. Jolicoeur, S. Ouvry, F. David (Springer- Verlag, 
Berlin and Les Editions de Physique, Les Ulis, 2000). 
H.A. Fertig, Phys. Rev. B. 40, 1087 (1989). 
E. Tutuc, M. Shayegan and D.A. Huse, Phys. Rev. Lett. 
93, 036802 (2004). 

L. Tiemann, J. G.S. Lok, W. Dietsche, K. von Klitzing, K. 
Muraki, D. Schuh, and W. Wegscheider, Phys. Rev. B 77 
033306 (2008). 

L. Tiemann, W. Dietsche, M. Hauser, and K. von Klitzing, 
New J. Phys. 10 045018 (2008); L. Tiemann, Y. Yoon, W. 
Dietsche, K. von Klitzing, and W. Wegscheider, Phys. Rev. 
B 80, 165120 (2009). 

K. Yang, K. Moon, L Zheng, A.H. MacDonald, S.M. 
Girvin, D. Yoshioka, and S.C. Zhang, Phys. Rev. Lett. 72, 
732 (1994) 

P. Bak, Rep. Prog. Phys. 45 587 (1982). 

C.B. Hanna, A.H. MacDonald and S.M. Girvin, Phys. Rev. 

B 63, 125305 (2001). 

G. D. Mahan, Many-Particle Physics (Plenum, New York, 
1990). 

Y. Yu and M. Cardona, Fundaments of Semiconductors: 
Physics and Matherial Properties (Springer 2001). 
P.R. Eastham, N.R. Cooper and D.K.K. Lee, Phys. Rev. 
B 80, 045302 (2009). 



